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The branched polymer growth model (BPGM) has been employed to study the kinetic growth of 
ramified polymers in the presence of impurities. In this article, the BPGM is revisited on the square 
lattice and a subtle modification in its dynamics is proposed in order to adapt it to a scenario closer 
to reality and experimentation. This new version of the model is denominated the adapted branched 
polymer growth model (ABPGM). It is shown that the ABPGM preserves the functionalities of 
the monomers and so recovers the branching probability b as an input parameter which effectively 
controls the relative incidence of bifurcations. The critical locus separating infinite from finite 
growth regimes of the ABPGM is obtained in the (b, c) space (where c is the impurity concen- 
tration). Unlike the original model, the phase diagram of the ABPGM exhibits a peculiar reentrance. 
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I. INTRODUCTION 

The statistical mechanics of polymers has become a 
very prolific research area. During the last decades, many 
models have been proposed to describe the configura- 
tional properties of linear and branched polymers and a 
variety of methods have been employed to study these 
systems. The classic model for a linear polymer in di- 
lute solution for at high temperature) is the self- avoiding 
walk (SAW) Q,0. On the other hand, connected clus- 
ters (lattice animals) provided a good model for branched 
polymers in the dilute limit [3( . 

The kinetic growth walk (KGW) was proposed as an 
alternative model to describe the irreversible growth of 
linear polymers [4-6]. Like the SAW, every KGW chain 
does not intercept itself. But whereas in a SAW the 
next step is randomly chosen from among all nearest- 
neighbor sites (excluding the previous one), in a KGW 
the choice is among the unvisited sites Q. Therefore 
the KGW is less sensitive to attrition. Besides, although 
both models span the same set of configurations, the 
N-step SAW chains are all equally weighted while the 
statistical weights of the KGW chains can be different. 
Nevertheless, it was shown that both models present the 
same critical exponents p. 

Later, Lucena et al. Q generalized the KGW in or- 
der to allow for branching as well as for impurities. This 
generalized model (which became known as the branched 
polymer growth model - BPGM) was found to exhibit an 
interesting phase transition (due to competition between 
hindrances and branching) separating infinite from finite 
growth regimes [jj. In the following years, several au- 
thors have studied the BPGM [8-15]. The topological 
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and dynamical aspects of the model were investigated 
Besides, the system was shown to achieve the criti- 
cality through a self-organization growth mechanism 
and to exhibit a transition from rough to faceted bound- 
aries for large values of the branching probability || El • 
The model was also studied through an exact enumera- 
tion of bond trees and ergodicity violation was discussed 
|l4j . The question of the universality class of the BPGM 
was also investigated. In contrast to the common belief 
that branched polymers belong to the universality class of 
lattice animals 0, , the study of the growth process 
in chemical /—space with estimates of structural expo- 
nents led to the proposal that the BPGM belongs to the 
universality class of percolation 0, . On the Bethe 
lattice, this proposal is based on analytical results 11]. 
A further analysis using finite-size scaling techniques led 
to the conclusion that the BPGM is not in the same uni- 
versality class of percolation in two dimensions [TEj . 

In the present paper, we revisit the BPGM on the 
square lattice. We point out that, although the input 
parameter b is named branching probability, it does not 
have effective control of the relative incidence of bifurca- 
tions (branches) in the simulated polymers. This seems 
to be an undesirable aspect of the model since the de- 
gree of branching is an important quantity that is usually 
measured in real ramified polymers. Besides, according 
to the growth rules of the BPGM, a monomer which is 
chosen to bifurcate (with probability b) but has only one 
empty nearest neighbor site is compelled to grow linearly 
thus changing its functionality from 3 to 2. In real experi- 
ments this change is not allowed since the concentrations 
of polymers with different functionalities are fixed. In 
order to preserve the functionalities of the monomers as 
well as to recover the meaning of the branching probabil- 
ity 6, we propose a subtle modification in the dynamics 
of the BPGM. So this new version of the model aims 
to adapt it to a more realistic scenario and is called the 
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adapted branched polymer growth model (ABPGM). In 
the following sections, we define the ABPGM and present 
the phase diagram in the (6, c) space. This diagram is 
found to exhibit a peculiar reentrance. We introduce the 
concept of frustration and compare the effective branch- 
ing rate b e f with the input parameter b for both the 
original and adapted models. Our results concerning the 
ABPGM are shown to be much closer to the ideal behav- 
ior b e f = b. Finally, we present a discussion based on the 
clusters topologies of both models. 



II. THE ADAPTED BRANCHED POLYMER 
GROWTH MODEL 

First of all, let us review the BPGM. Consider aLxL 
square lattice with a certain concentration c of sites ran- 
domly filled by impurities. At the initial time t = 0, the 
polymer starts growing from a monomer seed located at 
the center of the lattice towards a random empty nearest- 
neighbor site. At time t = 1, this chosen site is occupied 
by a monomer (growing tip) which now may bifurcate 
or follow in one direction (linear growth). The growth 
directions are randomly chosen among the available ones 
(those which lead to empty neighbors) in a clockwise way. 
At every time t ( t > 2), the process is repeated for all 
actual monomeric ends following the sequence of their 
appearances. If a polymer end has no empty nearest- 
neighbor site then it is trapped in a "cul de sac" and 
stops growing (it is then called a dead end). If only one 
adjacent site is available then the linear growth is oblig- 
atory. If at least two adjacent sites are empty then the 
growing end bifurcates with probability b or follows lin- 
early with probability 1 — b. In a particular experiment of 
the BPGM, the polymer growth is simulated according 
to the above rules. Depending on the values of parame- 
ters b and c, the polymer can either grow indefinitely or 
stop growing (die) at a finite time if all its current tips 
are dead ends. The experiment is finished either when 
the polymer touches the frontier of the lattice (infinite 
polymer) or when it dies (finite polymer). The ensem- 
ble over which averages are performed is constituted by 
a great number N e of experiments. 

Each polymer configuration of the BPGM can be iden- 
tified with a self-avoiding loopless graph (bond tree) al- 
though the reciprocal is not always true Thus the 
BPGM can be mapped into a particular subset of the en- 
semble of all possible bond trees. These tree graphs have 
been applied in models of branched polymers and can be 
classified by the number N of bonds and the number Nk 
of vertices with k bonds pi| . In a typical finite poly- 
mer configuration of the BPGM, one may find twofold 
and threefold sites representing monomers of function- 
alities 2 and 3, respectively. The monomer seed and 
dead ends represent monofunctional monomers. Tetra- 
functional units are impossible in this model once trifur- 
cations are not allowed. For a tree graph generated by 
the BPGM one has the following constraints |l8j : 
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^7V fc = 7V + i (1) 

k=l 

and 

3 

Y,kN k = 2N. (2) 

k=l 

The occurrence of bifurcations in the model deter- 
mines the relative incidence of trifunctional monomers 
(branches). This is basically the degree of branching 
of the polymer, an important quantity that is usually 
measured in real ramified polymers. However, although 
the parameter b is named branching probability, we shall 
see that it does not have effective control of the relative 
amount of branches (bifurcations) in the BPGM. Indeed, 
according to the rules of the BPGM, the parameter b is 
not the probability that any tip bifurcates but instead 
it is the conditional probability that a tip with two or 
more empty nearest neighbors bifurcates. Every tip with 
just one vacant nearest-neighbor site grows linearly with 
probability equal to one. From another viewpoint, it 
would also be desirable to consider that any monomer 
that is linking to the polymer has an effective probabil- 
ity b of being a trifunctional unit. Even if this were con- 
sidered in the BPGM, a problem would become evident. 
Any growing end chosen (with probability b) as a trifunc- 
tional monomer would not be able to bifurcate if only one 
vacant site were available. In this case, the functionality 
of the monomer would not be respected since (according 
to the BPGM) the growing end would forcibly follow a 
linear growth like a bifunctional monomer. 

In order to preserve the functionalities of the 
monomers we propose a modification in the dynamics of 
the BPGM. Besides to make the model a little more real- 
istic, our proposal turns the branching probability b into 
an effective control parameter of the relative frequency of 
bifurcations. So it adapts the model to a scenario closer 
to reality and experimentation. We denominate it the 
adapted branched polymer growth model (ABPGM). 

The ABPGM is defined just like the BPGM except for 
the following differences: 

1. In the process of polymerization, let a free end 
(a growing tip with at least one empty nearest- 
neighbor site) be a trifunctional monomer with 
probability b or a bifunctional monomer with prob- 
ability 1 — b. If the free end is a trifunctional 
monomer but there is only one empty nearest 
neighbor site then it stops growing and becomes 
a dead end. 

2. At every time unit, all current sites on the front 
of growth of the polymer are visited in a random 
sequence. 
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We remember that in the BPGM, a free end with just 
one empty nearest neighbor grows linear with probabil- 
ity one while in the ABPGM, since bifurcation is impos- 
sible, the free end stops and becomes a frustrated dead 
end. This frustration seems to be preferable and more re- 
alistic than forcing a trifunctional monomer to turn into 
bifunctional monomer as it occurs in the BPGM. So, in 
our proposal, the relative incidence of bifunctional and 
trifunctional monomers does not depend on the topol- 
ogy of the cluster anymore and is only controlled by the 
branching probability parameter. The second difference 
is also important: in the BPGM all tips are sequentially 
visited in a clockwise manner following the sequence of 
births whereas in the ABPGM they are visited in a ran- 
dom way. The sequential update of the growth front of 
the polymer (in the BPGM) simulates the formation of 
parallel chains in the infinite phase (for high b and low c) 
like in a crystallization process. This mechanism leads to 
a faceted-to-rough transition that has been already stud- 
ied Figure la shows a typical BPGM configura- 
tion with faceted boundary generated with parameters 
b = 0.5 and c = (for L = 51); the parallel ordering of 
chains is caused by the clockwise update of the growing 
tips. If we simulate the BPGM using the same set of pa- 
rameters but with a random update instead of the clock- 
wise update we get Figure lb. In this case, two subse- 
quent growing tips are probably located far apart so that 
they cannot produce parallel chains. The boundary is less 
faceted than before. Recently it has been shown that any 
deterministic growth order is non-ergodic in the sense 
that it spans only a subset of the space of all possible 
configurations [TJ- Moreover, for the present purposes, 
this clockwise update mechanism is undesirable since the 
development of parallel chains corresponds to an increas- 
ing incidence of bifunctional monomers due mainly to 
geometrical effects rather than to the probability 1 — b 
itself. The relative incidence of bifurcations (which will 
be defined as the effective branching rate b e f in section 
4) is much smaller than input parameter b. Indeed, we 
anticipate that b e f « 0.15 for the cluster of Figure la 
whereas b e f « 0.30 for the one of Figure lb; so the last 
rate (corresponding to a BPGM cluster generated with 
random update) is closer to the input value b = 0.50. 
The effective branching rate gets still closer to the input 
value if we simulate the ABPGM. Figure lc is a typical 
ABPGM configuration (generated with those same pa- 
rameters b = 0.5 and c = 0) and presents b e f ^0.43; the 
faceted front of growth seems to disappear and vacancies 
can be found inside the cluster. 

The main feature of the BPGM is the phase transi- 
tion separating finite from infinite growth regime. Lu- 
cena et al. Q have defined a critical branching proba- 
bility b c where the mean size of finite polymers diverges 
as L — > oo. The critical value depends on the impurity 
concentration c and separates the finite phase (b < b c (c)) 
from the infinite one (b > b c (c)). The probability that 
a polymer grows indefinitely is null for small values of b 
and increases abruptly in the region of b « b c (c). This 
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FIG. 1: Typical configurations simulated on a square lattice 
of size L = 51 with parameters b = 0.5 and c = according 
to the BPGM rules with clockwise update of growing tips 
(a), BPGM with random update (b) and ABPGM (c). The 
corresponding effective branching rates are b e f ~ 0.15, b e f ~ 
0.30 and b e f ~ 0.43, respectively. 



probability was estimated through the fraction Poo of in- 
finite polymers in the ensemble of configurations of the 
BPGM 71. In Fi gure 2a, we reproduce typical plots of 
Poo versus b at different values of c (with N e = 10 2 ex- 
periments and L = 1501) for the BPGM. In Figure 2b, 
we present plots of Poo versus b at some values of c cor- 
responding to simulations of the ABPGM with N e = 10 4 
experiments and L = 1501. For c = 0, the behavior of 
Poo (b) is analogous to that of the BPGM but now the 
threshold is little bit higher (b « 0.06 for the ABPGM 
while b « 0.055 for the BPGM). This difference increases 
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with c. The most interesting characteristic of the present 
model may be observed when c = 0.18. For this value, the 
curve raises at b « 0.32, then presents a plateau (where 
P<x>(b) « 0.7) and finally falls! By a finite-size scaling 
analysis we have verified that the height of the plateau 
does not change significatively in the limit L — > oo. For 
c = 0.19, the behavior of Poo(b) is gaussian shaped. This 
means that, for certain impurity concentrations, as b in- 
creases the system goes from a finite to an infinite phase 
and then becomes finite again! Indeed, this reentrance 
is confirmed in the next section when we determine the 
phase diagram of the ABPGM through an analysis of the 
correlation length. 
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FIG. 2: Plots of the fraction of infinite polymers 

Poo versus b 

for the BPGM (a) and ABPGM (b). 



III. THE PHASE DIAGRAM 

The mean size of finite polymers is a measure of the 
correlation length £ of the system. If a finite polymer 
is generated during a simulation, the sizes l x and l y of 
the smallest rectangle containing the cluster can be de- 
termined. The correlation length can then be calculated 

1/2 

as £ = ((Wy) ) where the average is performed over 
all experiments with finite polymers Q. We show typ- 
ical plots of £ versus b corresponding to simulations of 
the BPGM (in Figure 3a) and ABPGM (in Figure 3b) 



for different values of c and size L = 1501 (with 10 2 and 
10 4 experiments respectively). Each plot of £ versus b ex- 
hibits only one maximum except for plots of the ABPGM 
with c > 0.175 where two peaks are detected! It can be 
verified that all peaks of £ do diverge when L — > oo. 
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FIG. 3: Plots of the correlation length £ of finite polymers 
versus the branching probability b for several impurity con- 
centrations c in the cases: (a) BPGM and (b) ABPGM. 

Let us first explain the behavior of £ for the BPGM. 
For a fixed value of c and as b increases, the polymer 
is more likely to escape from steric hindrances and im- 
purities so that the mean size grows to its highest value 
at some b = bo(L). Above this point, the system is de- 
fined to be in the infinite growth regime (of course, the 
true critical point is obtained in the thermodynamic limit 
b c = liniL^oo bo(L)). As b continues to increase, the frac- 
tion of infinite polymers grows and the finite polymers get 
smaller so that £ decreases. As higher impurity concen- 
trations hindrance the growth, bo increases with c. The 
critical line of the BPGM is the locus (b c ,c) on which £ 
diverges Q and is shown in Figure 4 (dashed line). 

Regarding the ABPGM, the same reasoning can ex- 
plain the maximum of £ (or the first maximum when 
there are two peaks). But now this peak is located on 
a higher branching probability b\ > bo that compensates 
the occurrence of frustrated dead ends. Just above bi 
the system enters the infinite growth regime where it re- 
mains unless a second peak appears (at b = 62 > b\ for c 
> 0.175). In the latter case the system returns to the fi- 
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nite phase! This reentrance from infinite to finite growth 
regime is a peculiar feature of the ABPGM. For this mod- 
ified model, there is a certain range of values of c where 
it is very probable that all free ends become frustrated 
(and stop growing) for a sufficiently large b so that the 
polymer growth is finite again. The ABPGM phase dia- 
gram is also shown in the Figure 4. The reentrant phase 
only exists for c in the small interval [0.175,0.195]. 
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FIG. 4: The phase diagram of the models BPGM (dashed 
line) and ABPGM (solid line). Each critical line separates 
the infinite phase (at low values of c) from the finite one (at 
higher values of c). 



critical line, P(N) decays exponentially with TV. Lucena 
et al. 7] have found a similar behavior in the BPGM. 

IV. THE FRUSTRATION AND THE 
EFFECTIVE BRANCHING RATE 

The frustration is a desirable event which distinguishes 
the ABPGM from the original model. It is defined as the 
interruption of the growth of any free end which was cho- 
sen (with probability b) as a trifunctional monomer but 
is unable to bifurcate since only one nearest neighbor site 
is available. It prevents such a free end to continue its 
growth linearly like a bifunctional monomer (as occurs in 
the BPGM) and consequently controls the relative inci- 
dence of branches. 
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FIG. 5: Log- log plot of the distribution P(N) of finite poly- 
mers with N bonds for simulations of the ABPGM (with 
JV" e = 10 4 and L = 1501) at some points (6, c) on the crit- 
ical line. 

We have also measured the fraction P(N) of finite 
polymers with N bonds. In Figure 5, we have a log- 
log plot of the polydispersion distribution P(N) of the 
ABPGM on the critical line (running 10 4 experiments in 
a L = 1501 lattice). The three sets of data correspond 
to the critical points: b\ « 0.06 and c = (black dots); 
bi « 0.32 and c = 0.18 (open circles); b^ ~ 0.93 and 
c = 0.18 (squares). The data are fitted by a straight line 
with slope —0.75. So, on the critical line, P(N) decays 
with N as a power law. We have verified that outside the 



FIG. 6: The frustration rate f versus b for increasing values 
of c in the ABPGM (L = 1501 and iV e ~ 10 3 ). 

Clearly, each polymer configuration of the ABPGM is 
also a bond tree which can be classified by the numbers N 
(of bonds) and Nk (of k— functional units). Besides, every 
dead end is a site which ceased to grow and so represents 
a monofunctional vertex; it can be subclassified as either 
a trapped site (if it is in a "cul de sac") or a frustrated 
dead end (otherwise). We denote by Ni the number of 
frustrated dead ends. For the BPGM, Ni = since any 
free end never stops. In order to measure the relative 
incidence of frustrated ends among those free sites chosen 
as trifunctional monomers, we define the frustration rate 
as 

f=(^-) (3) 

\n 1 + n 3 / 

where the average is performed over all experiments. For 
infinite polymers, the current sites on the front of growth 
should not be considered in the computation since their 
functionalities are undetermined. 

According to the prior definition, / = for any values 
of b and c in the BPGM. We remark that the null frustra- 
tion of the BPGM does not mean that all free sites chosen 
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to bifurcate succeed but only that when they fail they are 
transformed in monomers with functionality two. 

For the ABPGM, the frustration rate increases with 
b and c as it is shown in Figure 6. Indeed, both the 
excluded volume due to self- avoidance (which increases 
with b) and the impurities diminish the chance of suc- 
cess of any free end, so / is a monotonically increasing 
function of b and c. 
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FIG. 7: The behavior of the effective branching rate b e f versus 
the input parameter b for simulations of the BPGM (a) and 
ABPGM (b) on a large L = 1501 square lattice (with 10 2 and 
10 4 experiments respectively). The dashed lines represent the 
ideal behavior b e f = b. 

The effectiveness of the input parameter b is evaluated 
by comparing it with the relative frequency of branches 
in a polymer configuration. For this purpose, we define 
the effective branching rate b e f as the ratio between A3 
and N2 + A3 averaged over all experiments: 



N 3 



N2 + A3 



(4) 



The plots of b e f versus b at different values of c for sim- 
ulations of the BPGM and ABPGM (on a large L = 1501 
square lattice) are shown in Figures 7a and 7b, respec- 
tively. The ideal behavior b e f = b is indicated as the 
dashed lines. Regarding the BPGM, there is evidently 



a large discrepancy between the input parameter b and 
the output rate b e f. For c = 0, b e f increases with b up 
to a maximum (at b ~ 0.2) and then decreases to zero as 
b —> 1. This behavior is explained as follows: for small 6, 
self- avoidance is reduced so that most sites trying to bi- 
furcate succeed; but as b increases, parallel linear chains 
(like those of Figure la) are forcibly generated due to 
both the increasing excluded volume and the clockwise 
update. For higher values of c, the presence of impuri- 
ties hampers the formation of linear chains and thus b e f 
decreases slower (as a concave function). Anyway, for 
the BPGM, the discrepancy between b and b e f gets more 
pronounced as b — » 1. On the contrary, for the ABPGM, 
b e f always increases with b as a convex function (for any 
value of c) and the difference b — b e f is very small. In fact, 
for any c, the ratio b e f/b is approximately equal to 1 for 
small 6, decreases until about 0.8 for intermediate b and 
then returns to 1 as b —> 1. Such results corroborate the 
assertion that the input parameter b controls the relative 
incidence of branches in the adapted model. 
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FIG. 8: Typical graphs of both the BPGM (a) and ABPGM 
(b) for b = 0.2 and fixed c = 0.18. 



V. CONCLUSIONS 

The branched polymer growth model (BPGM) was orig- 
inally proposed | as a generalization of the kinetic 
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b e f is very much closer from the input parameter b in the 




growth walk Q in order to include the possibility of ram- 
ification of the polymer as well as the presence of impu- 
rities in the medium. The model was found to exhibit 
a finite-infinite transition due to competition between 
branching and hindrances. 

In this paper, we have proposed an alteration in the 
dynamics of the BPGM so as to adapt the model to 
an experimental realism. We have called it the adapted 
branched polymer growth model (ABPGM). The main 
difference between our proposal and the original model 
regards to the growth mechanism of a monomer which 
is chosen to bifurcate (with probability b) but has just 
one empty nearest neighbor site. In the BPGM, such a 
monomer is transformed into a bifunctional unit so that 
it grows linearly (with probability one). In our adapted 
model, that monomer stops and becomes a frustrated 
dead end. This frustration reveals as preferable and more 
realistic than changing the monomer functionality from 
3 to 2 (as it occurs in the BPGM). This subtle change 
in the algorithm together with a random update of the 
growing ends lead to the formation of polymers with new 
topological patterns and adjusted degrees of branching. 
Indeed, we have shown that the effective branching rate 



FIG. 10: Typical graphs of both the BPGM (a) and ABPGM 
(b) for b = 1 and fixed c = 0.18. 



trant phase in the small interval 0.175 < c < 0.195. At 
this instance, we compare some typical graphs of both the 
BPGM and ABPGM at the fixed impurity concentration 
c = 0.18 for three increasing values of b: 0.2, 0.6 and 1 
(Figures 8, 9 and 10, respectively). For the BPGM, all 
the three corresponding configurations are infinite clus- 
ters whose boundaries change from rough to faceted as b 
increases. On the other hand, the typical ABPGM graph 
for b = 0.2 is a finite cluster (due to the occurrence of 
both trapped sites and frustrated dead ends); if b = 0.6 
the cluster is infinite (since here higher branching over- 
comes dead ends) and if b = 1 the cluster is finite again 
(due to a high frustration rate). 

We hope that the ideas presented here can add a new 
degree of physical understanding to the model as well as 
they can approximate it to the real ramified polymers. 

This work has been supported by Brazilian Agen- 
cies CAPES, CNPq and FAPESP (through grant No. 
96/05387-3). 
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